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The implications of the first AMS-02 p/p data for the propagation of cosmic rays and 
the properties of dark matter (DM) are discussed. Using various diffusive re-acceleration 
(DR) propagation models, one can derive very conservative upper limits on the DM 
annihilation cross sections. The limits turned out to be compatible with that from the 
Ferm-LAT gamma-ray data on the dwarf spheroidal satellite galaxies. The flattening of 
the p/p spectrum above ~ 100 GeV in the current data still leaves some room for TeV 
scale DM particles. More antiproton data at high kinetic energies are needed to constrain 
the properties of the DM particles. 
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Cosmic-ray antiparticles, such as positrons and antiprotons play important roles 
in the indirect search for dark matter (DM) in the Galactic halo. The Alpha Mag¬ 
netic Spectrometer (AMS-02) is measuring such cosmic-ray charged particles with 
unprecedented accuracies. So far the anomalous rise in the positron fraction pre¬ 
vious reported by PAMELAPI^ and Fermi-LATl^ has been confirmed by AMS-02 
with higher accuracy and extended to higher energies,® which has triggered exten¬ 
sive theoretical studies on possible explanations including halo DM annihilation or 
decay (for recent global analyses on AMS-02 data, see e.g. Ref^^^^ ). Antiprotons 
are highly expected from DM annihilation in many DM models, which is unlikely 
to be generated from the nearby pulsars. 

Recently, the AMS-02 collaboration has released the first preliminary result of 
the cosmic-ray antiproton to proton flux ratio p/p The measured kinetic energies 
of the antiprotons have been extended to ^ 450 GeV. Although the spectrum of p/p 
at high energies above 100 GeV tend to be relatively flat, within uncertainties the 
AMS-02 data are consistent with the background of secondary antiprotons, which 
can be used to set stringent upper limits on the dark matter (DM) annihilation 
cross sections, especially for high mass DM particles. The constraints on the DM 
properties from antiprotons have been investigated previously before AMS-02 ( see 
0 gjISllI] y jjj talk, we briefly summarise our work on the implications of the 
new AMS-02 p/p data for constraining the annihilation cross sections of the DM 
particles in various propagation models and DM profiles. The details of the analysis 
can be found in Ref.^^ 

In the diffusion models of cosmic-ray propagation, the Galactic halo within which 
the diffusion processes occur is parametrized by a cylinder with radius Rh = 20 — 30 
kpc and half-height = 1 — 20 kpc. The diffusion equation for the cosmic-ray 
charged particles reads 

pV:-|(v-y,)v^' 

- —i’- —'4’ + q{r,p), (1) 

Tf Tr 

where ip(r,p,t) is the number density per unit of total particle momentum. For 
steady-state diffusion, it is assumed that dip/dt = 0. The number densities of 
cosmic-ray particles are assumed to be vanishing at the boundary of the halo. 
The energy dependent spatial diffusion coefficient D^x is parametrized as Dxx = 
I3Dq {p/poY ; where p is the rigidity of the cosmic-ray particle. The power spectral 
index 6 can have different values 5 = ^ 1 ( 2 ) for p below (above) a reference rigidity 
Po- Dq is a normalization constant. The convection term in the diffusion equation 
is related to the drift of cosmic-ray particles from the Galactic disc due to the 
Galactic wind. The diffusion in momentum space is described by the reacceleration 
parameter Dpp which is related to the Alfven speed Va of disturbances in the hy- 
drodynamical plasma.!^ The momentum loss rate is denoted by p, and Tf{Tr) is the 
time scale for fragmentation (radioactive decay) of the cosmic-ray nuclei. 
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The spectrum of a primary source term for a cosmic-ray nucleus A is assumed to 
have a broken power low behaviour dqA(j>)/dp oc {p/pAsV^ with 7 ^ = 7ai( 7A2) for 
the nucleus rigidity p below (above) a reference rigidity pAs ■ The spatial distribution 
of the primary sources is assumed to follow that of the pulsars and is taken from 
Ref.^ The background antiproton is assumed to only have the secondary origin, 
namely, they are created dominantly from inelastic pp- and pA-collisions with the 
interstellar gas. The corresponding source term reads 

gsec(p) = Pen, y] f dp'Ahflnp{p), (2) 

where Ui is the number density of the interstellar hydrogen (helium), np is the num¬ 
ber density of primary cosmic-ray proton per total momentum, and dai{p,p')/dp' 
is the differential cross section for p + H(He) p + X. In calculating the antipro¬ 
tons, inelastic scattering to produce “tertiary” antiprotons should be taken into 
account. The primary source from the annihilation of Majorana DM particles has 
the following form 

= . ( 3 ) 

where (av) is the velocity-averaged DM annihilation cross section multiplied by DM 
relative velocity. p{r) is the DM energy density distribution function, and dN^^'>/dp 
is the injection energy spectrum of antiprotons from DM annihilating into SM final 
states through possible intermediate states X with px the corresponding branching 
fractions. The interstellar flux of the cosmic-ray particle is related to its density 
function as $ = vip{r,p)/{4:Tr). At the top of the atmosphere (TOA) of the Earth, 
the fluxes of cosmic-rays are affected by solar winds and the helioshperic magnetic 
field. This effect is taken into account using the force-field approximation which 
involves the Fisk potential We shall take p = 550 MV in numerical analysis. 

We solve the diffusion equation of Eq. o using the publicly available code GAL- 
PROP v 542 SH 20 | which utilizes realistic astronomical information on the distribution 
of interstellar gas and other data as input, and considers various kinds of observ¬ 
ables in a self-consistent way. We start with the so-called “conventional” diffusive 
re-acceleration (DR) modeP^^ which is commonly adopted by the current exper¬ 
imental collaborations as a benchmark model for the astrophysical backgrounds. It 
is useful to consider this model as a reference model to understand how the DM 
properties could be constrained by the AMS-02 data. Then we consider three repre¬ 
sentative propagation models selected from a large sample of models obtained from 
a global Bayesian MCMC fit to the preliminary AMS-02 proton and B/C data using 
the GALPROP code.^They are selected to represent the typically minimal (MIN), 
median (MED) and maximal (MAX) antiproton fluxes within 95% GL, correspond¬ 
ing to the region enveloping 95% of the MGMG samples with highest likelihoods in a 
six-dimensional parameter space. Note that the GALPROP based “MIN”, “MED” 
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FIG. 1. Predictions for the p/p ratio from the four propagation models. The data from AMS-OJ^ 
and PAMELAP^are shown. See text for detailed discription. 


and “MAX” models used in this work are different from and complementary to that 
given in RefP 

The predictions for the background of the p/p flux ratio in these models are 
shown in Fig. [TJ The “MIN”, “MED” and “MAX” models are highly degenerate in 
the background p/p ratio. Compared with these models, the “conventional” model 
predicts more low energy antiprotons but at high energies above ~ 500 GeV, the 
predicted antiprotons are much less. In all the four DR propagation models, below ~ 
10 GeV the GALPROP based calculations underpredict the p/p flux ratio by ^ 40%, 
which is a known issue. The agreement with the low energy p data can be improved 
by introducing breaks in diffusion coefficients,^^ “fresh” nuclei componenlP^or a DM 
contribution.^^ The predictions for low energy p/p ratio can also be easily modified 
by introducing an independent Fisk potential (j) for p and an energy-dependent 
overall normalization factor as discussed in Ref.^^ We instead use these DR models 
to derived very conservative upper limits on the annihilation cross sections of light 
DM particles. Note, however, that in the DR propagation models, the background 
predictions agree with the AMS-02 data well at higher energies ^ 10 — 100 GeV, 
which can be turned into stringent constraints on the nature of heavy DM particles. 

We consider three reference DM annihilation channels XX where XX = 

gg, bb and W~^W~. The energy spectra of these channels are similar at high energies. 
The main difference is in the average number of total antiprotons Xx per DM 
annihilation of each channel. The injection spectra /dp from DM annihilation 

are calculated using the numerical package PYTHIA v8. 175.1^ 

We shall first derive upper limits on DM annihilation cross section as a function 
of DM particle mass, using the frequentist x^-analyses. All of the 30 data points 
of the AMS-02 p/p data are included in calculating the limits. In Fig. [21 we show 
the obtained upper limits on the cross sections for DM particle annihilation into 
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FIG. 2. Upper limits on the cross sections for DM particle annihilation into bb final states from 
the AMS-02 p/p data in the “conventional” (upper left), “MED” (upper right), “MIN” (lower left) 
and “MAX” (lower right) propagation models. Four DM profiles NFWp^Isotherma®! , EinastJlSl 
and MoorJSIlll are considered. The upper limits from the Fermi-LAT 6-year gamma-ray data of 
the dwarf spheroidal satellite galaxies of the Milky Way are also shownBS] The horizontal line 
indicates the typical thermal annihilation cross section (cru) = 3 X 10“^®cm®s“^. 


bb final states from the AMS-02 p/p data in the “conventional”, “MED”, “MIN” 
and “MAX” propagation models. Eour different DM profiles: NEW,I^ Isothermal,!^ 
EinastcP^ and Moor^SHH are considered. As can be seen, the upper limits as a 
function of show some smooth structure for all the final states and DM profiles. 
The limits tend to be relatively stronger at « 300 GeV, which is related to 
the fact that the background predictions agree with the data well at the antiproton 
energy range ^ 20 — 100 GeV. Eor a comparison, the upper limits from the Fermi- 
LAT 6-year gamma-ray data of the dwarf spheroidal satellite galaxies of the Milky 
Waj®! are also shown in Fig. [2] In the “conventional” model, the upper limits from 
the AMS-02 pjp data are found to be compatible with that derived from the Fermi- 
LAT gamma-ray data for rriy^^ > 300 GeV. This observation holds for most of the DM 
profiles. In the “MED” model, the constraints are relatively weaker, which is related 
to the under prediction of low energy antiprotons in this model and the limits are 
more conservative. Eor an estimation of the uncertainties due to the propagation 
models, from the “MIN” model to the “MAX” model, we find that the variation of 
the upper limits is within about a factor of five. 

For the W~^W~ final states, the resulting limits are shown in Fig. [3l In the 
“conventional” propagation model, the constraints from AMS-02 p/p data turn out 
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FIG. 3. The same as Fig. (2] but for DM annihilation into W'^W final states. 






FIG. 4. The same as Fig. (2] but for DM annihilation into qq final states. 
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to be more stringent than that from the Fermi-LAT gamma-ray data for all the four 
DM profiles when the DM particle mass is above ^ 300 GeV. Again we find that 
the variation of the upper limits from the “MIN” to the “MAX” model is within 
a factor of five. The result for the qq final states is shown in Fig. |4l Compared 
with the case of W~^W~ and bb, the constraints on the qq final states are the most 
stringent. For all the three final states, we find that the allowed DM annihilation 
cross section is below the typical thermal cross section for < 300 GeV in the 
conventional propagation irrodel with Einasto profile. 




FIG. 5. Left) Spectrum of p/p flux ratio from DM annihilating into bb flnal states with = 
58.5 GeV and (av) = 2.16 X 10“^® cm®s“^ obtained from a fit to the whole AMS-02 p/p data.ESl 
The “conventional” background model and the Einasto DM profile are assumed. Right) The same 
as left, but for the fit with qq final state with the best-fit values = 35 GeV and (av) = 
0.86 X 10-2® cm®s-L 


As shown in Fig.[Tl compared with the AMS-02 data the GALPROP DR models 
predict fewer antiprotons at low (< 10 GeV) and very high (> 100 GeV ) ener¬ 
gies. Without a robust estimation of the theoretical uncertainties, it is too early 
to claim any excesses in the p/p data. We nevertheless consider what would be 
the implications for DM if such a trend in the observations is confirmed by future 
analyses. The low energy data would allow for a non-vanishing DM annihilation 
cross section. Eor instance, in the “conventional” propagation model, for =10.1, 
35.0 and 75.8 GeV, the best-fit values are {av) = 3.6 x 10“^^, 1.14 x 10“^®, and 
2.79 X 10“^® cm®s-^, respectively, if the DM prohle is Einasto, and the DM par¬ 
ticles annihilate dominantly into bb final states. If both and (av) are allowed 
to vary freely, the best-fit DM particle masses and annihilation cross sections are 

= 58.5 (35.0) GeV and (av) = 2.16 (0.86) x 10“^® cm^s”^ for DM annihilating 
into bb {qq) final states. In Eig. [5l we show the calculated spectra of p/p flux ra¬ 
tio from the best-fit DM particle masses and cross sections. The figure shows that 
the low energy p/p data are well reproduced by including such a DM contribution, 
except for the data point with kinetic energy below 1 GeV. 

As shown in Eig. [T] the spectrum of the AMS-02 p/p ratio tends to be flat to¬ 
ward high energies above ^ 100 GeV. This trend, if confirmed by the future AMS-02 
data, is not expected from the secondary production of antiprotons, and raises the 
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FIG. 6. (Upper left) values of Xmin ^ function of DM particle mass irix from a fit to t he AM S-02 
p/p data ( with kinetic energy above 20 GeV ) in the “conventional” propagation model with 

the DM profile fixed to EinastoBHl Three annihilation channels bb, qq and W'^W~ are considered. 
(Upper right) predicted p/p ratio in the case of background (“conventional” model) plus a DM 
contribution with = 6.5 TeV, (crv) = 1.9x 10~^‘^cm^s“^, and annihilation final states W^W~. 
The flux ratio of antiproton from DM to the proton from the background Pdm/pbg is shown as 
the dashed line. The data from AMS-Ofl^ and PAMELAEH are also shown. (Lower left) the same 
as the upper right, but for the bb channel with = 10.9 TeV and (av) = 3.4 X 10“^^ cm®s“^. 
(Lower right) the same as the upper right, but for the qq channel with = 10.9 TeV and 
(av) = 3.3 X lO-^'* cm^s-l. 


interesting question whether this would leave some room for a heavy DM contri¬ 
bution, similar to the case of the AMS-02 positron fraction.^l^lSHII] Xo explore this 
possibility, we perform an other fit using the p/p ratio data above 20 GeV (including 
15 data points in total) in order to avoid the theoretical uncertainties in the low 
energy region. The obtained Xmin a fnnction of for the bb, qq and W~^W~ 
final states in the “conventional” propagation model with Einasto DM profile are 
shown in Eig. [6l For the three final states the values of Xmin decrease almost mono- 
tonically from ^ 21 to ^ 5 with an increasing DM particles mass from 100 GeV 
to 10 TeV, but the x^-curves become gradually flat toward high DM masses. Only 
for the W~^W~ channel, there exists a shallow local minimal at around 6.5 TeV 
with low statistical significance. From the x^-curves, one can see that the DM par¬ 
ticles mass is restricted to be above ^ 2 TeV at 2a. For an illustration pnrpose, 
we show in Fig. [6] the predictions for the p/p ratio in the “conventional” back¬ 
ground model with a DM contribution. The DM particles masses and annihilation 
cross sections chosen to be m,, = 6.5 TeV, {av) = 1.9 x 10“^^cm^s“^ for W^W~, 
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rriy^ = 10.9 TeV, {av) = 3.4 x 10“^'’’ cm^s“^ for bb channel, and = 10.9 TeV 
and {(Tv) = 3.3 x 10“^^ cm^s“^ for qq channel. Note that these values are not from 
the best-fit values. We conclude that introducing a DM contribution can improve 
the agreement with the AMS-02 p/p data with kinetic energy above 100 GeV, but 
the statistics is not high enough to determine the DM properties such as its mass 
and interaction strength. As can be see in Fig.[T] the possible “excess” is located at 
the kinetic energy range 100 — 450 GeV where the secondary backgrounds from the 
four propagation models are similar. However, beyond ~ 450 GeV, the p/p from 
the “conventional” model drops quicker than that in the other propagation models. 
The future high energy antiproton data will be very important not only in probing 
DM but also in constraining the background models. 

In conclusion, we have explored the implications of the first AMS-02 p/p data on 
constraining the annihilation cross sections of the DM particles in various propaga¬ 
tion models and DM profiles. We have derived the upper limits using the GALPROP 
code and shown that in the “conventional ” propagation model with Einasto DM 
profile, the constraints can be more stringent than that derived from the Eerm-LAT 
gamma-ray data on the dwarf spheroidal satellite galaxies. Making use of the typ¬ 
ical minimal, median and maximal models obtained from a previous global fit, we 
have shown that the uncertainties on the upper limits is around a factor of five. 
The future more precise AMS-02 data can help to reduce the uncertainties in the 
derived upper limits. The analysis in this work has some overlap with that irP^ . 
Note that although the conclusions are similar, the analysis in this work is based on 
the fully numerical GALPROP code, while that irP^is based on the two-zone diffu¬ 
sion model with (semi)-analytical approach. Similar discussions on the DM matter 
contributions can be found in Refs.^^^^ 

This work is supported in part by the National Basic Research Program of 
Ghina (973 Program) under Grants No. 2010GB833000; the National Nature Sci¬ 
ence Eoundation of Ghina (NSEC) under Grants No. 10905084, No. 11335012 and 
No. 11475237; The numerical calculations were done using the HPC Cluster of 
SKLTP/ITP-CAS. 
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